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1. Introduction 



To establish the formation of a thermalized quark-gluon plasma during a heavy-ion 
collision it is crucial to determine the time scales involved in the equilibration of 
' the initial fireball. The hydrodynamic description of the experiments suggests early 

thermalization [ ^ , which seems to contradict traditional perturbative estimates 
[ 12] . The resolution of this puzzle has spurred the development of techniques to 
understand the microscopic dynamical processes leading to equilibration. 

A powerful formalism to study out-of-equilibrium dynamics from first principles 
is the 2PI effective action. It provides an exact representation of a given theory 
in terms of a functional depending solely on the connected one- and two-point 
functions. Approximate equations describing the time evolution of these correlators 
can be derived from a variational principle on the functional. The main interest 
of the 2PI effective action as a framework to obtain evolution equations is the fact 
that the global symmetries of the theory are respected [|n|. In particular, the 
evolution equations guarantee energy conservation, which is an important feature 
when studying out-of-equilibrium processes. They also seem to be free of the secular 
instabilities that plague traditional perturbative approaches. 
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In recent years, approximations based on the 2PI effective action have been 
apphed succesfuUy to the study of far-from-equilibrium dynamics and equihbration. 
Most studies have been carried out in scalar theories, both in 1+1 [01 |S1 IH] and 
2+1 dimensions [[7j, and in Yukawa theory in 3+1 dimensions [ 8 . The application 
of 2PI effective action methods to gauge theories, however, is not straightforward 
due to a residual dependence on the choice of gauge condition [ 9 and the difhculty 
of implementing properly a numerical discretisation. In this work we extend to 
3+1 dimensions previous studies of equilibration in scalar theory, investigating, in 
addition, the behavior in the broken phase. More information can be found in [ll(J|. 

2. 2PI Effective action and evolution equations 

For the scalar theory given by the action in the form 

^ Jd'x I d'y^ip{x)G,\x,yMy) - ^^(x)^ (1) 
with Gq^{x, y) ~ — m?)6{x^ y), the corresponding 2PI effective action is given 

by m 

r2Pi[0, G] - S[^] - ^TrlnG + ^T^t[{Go' - G'') ■ g] + $[0, G] (2) 

Here (j) f^nd G correspond, respectively, to the full connected one- and two-point cor- 
relation functions. The functional $ consists of an inifinite series of closed skeleton 
diagrams 



(3) 



where the lines correspond to G, the crosses to (p and the vertices to —iX. 

Evolution equations for the one- and two-point functions are obtained by a 
variational principle on the effective action r2Pi, conveniently formulated using a 
real-time contour for the time integrations. One obtains 



5r 



^-^^!^=O^Sc{x,y)= f d^zGo\x,z)G{z,y)+i f d^zJ:{x,z)G{z,y) (4) 



SG 

+ + (5, 

with C being the integration time contour (see for instance [.12,'). The self-energy 
y) is given by 

SG{y,x) 

The exact evolution of and G would be determined by solving these equations 
including the infinite series of diagrams in $[0, G]. This is not possible to do, so one 
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Truncation 



Hartree approximation j ~'~ ^ 3 

Two-loop approximation j + I Q + II " 

"Basketball" approximation j ^ ^ 3 T2 " " + Js 



Table 1. Truncations of the 2PI effective action. 

considers approximate evolution equations that result from taking only a finite set 
of diagrams, i.e. by truncating $[0, G] up to some order in a convenient expansion 
parameter. In our study we consider various truncations of the loop expansion of 
r2Pi, summarized in Tabled 

In a real-scalar theory, the two-point functions defined on the time contour C 
can be written in terms of two independent components, which we take as 



Gcix,y) = F{x,y) - -signc(a;o - yQ)p{x,y). 



(7) 



The functions F and p contain, respectively, statistical and spectral information 
about the system. They are real and satisfy the symmetry properties F{x,y) = 
F(y,x) and p{x,y) = —p{y,x), which make them very convenient for numerical 
implementation. In terms of these functions, the evolution equations look like 

r^o r ryo r 

[dl + M^{x)] F{x, y)^ dzo / d^z z)F{z, y) - dza / d^z ^^{x, z)p{y, z) 



(8) 



[dl + M\x)] p{x, y) = / dzo d?z ^P{x, z)p{z, y), 



Xq 



[dl + M^ix)] 4>{x) = -c^{xf + J dzo J d^z ^P{x, z)4>{z), 

where M'^{x) = m? + -I- ^F{x,x). For the Basketball approximation^ we 

have 



cj){x)(i)(y) 



F^{x,y) 



P^{x,y) 



\2 

+ —F(x,y) 
6 



F\x,y) 



3p2(x,y) 



T.P{x,y) = \^(l){x)(l){y)[F{x,y)p{x,y)\ + —p{x,y) 



iF\x,y) 



p'^{x,y) 



(9) 
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■p{x,z) 



\F{x^,f.P±A 



We see from (jS)) that the evolution oi F, p and (p is determined by their values at 
previous times, which enter the equation in the form of memory integrals. For the 
two-loop approximation, the memory kernels are given by the first term in the RHS 
of 0, while for the Hartree approximation they vanish alltogether. 
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To solve set of coupled evolution equations © we need to specify a set of initial 
conditions. We consider a spatially homogeneous situation so that the correlators 
can be written in terms of their mode functions F\^{t,t') and pk{t,t'). The initial 
conditions are then given by the values and derivatives of the mode functions F, p 
and (j) at initial time. In the case of the spectral function p, these are specified by 
the property 

Pk(i,i) = 0, dMi.t\^,, = l. (10) 
The initial conditions for F are taken of the form 



dtdt.F^,{t,t') 



1 



(11) 



which are identical to the equilibrium free theory case, but with some arbitrary 
distribution function rik- 

The evolution equations ^ are solved numerically by discretizing the theory 
on a space-time lattice with spacings at a. The lattice provides a cut-off and reg- 
ularizes the ultraviolet divergences present in the continuum limit, which are to be 
dealt with by renormalization. The continuum renormalization of approximations 
based on the 2PI effective action has been recently studied in detail (see for instance 
[ llHj V For our purpose it is sufficient to use an approximate renormalization that 
ensures that the relevant length scales in the simulations are larger than the lattice 
spacing a. For all the truncations of table we take a perturbative version of the 
renormalization needed just for the Hartree case. 



3. Symmetric Phase 

We consider first the evolution of the system in the symmetric phase. We initialize 
the system far from equilibrium with a distribution of the form 

nk = rye(kLx - k2)e(k2 - k^iJ, (12) 

which means that only modes with momenta in the range k^^jj^ < k^ < k^^^^ are 
occupied. In the following we shall consider three diferent distributions: Tl, T2 and 
T3, all with identical energy and Tl and T2 with similar initial particle number. 
The numerical evaluation of the evolution equations © is performed on a 16'^ spatial 
lattice with spacing am ~ 0.7 (in units of the renormalized mass m), time spacing 
atm = 0.07 and coupling A = 6.^ To monitor the evolution of the system towards 

^The memory integrals in jSJ arc cut to a maximum time extent mtcut ~ 28. Contributions to 
the kernels from previous times were checked to be negligible. 



Kinetic and Chemical Equilibration in Scalar (^^ Theory 



5 




Fig. 1. Evolution in time of the occupation numbers, Tik vs. ujy^ (left), and dis- 
persion relation vs. (right), for an initial condition with (77, k^^^, k^^^) = 
(2,2.04m2,6.12m2) (Tl). We display the results of the "Basketball" (black dots) 
and the Hartree approximation (green/grey dots). 

equilibrium we study the behavior of an effective quasiparticle distribution function 
and a dispersion relation, defined as 

n^{t) + c^^dtdvF^,{t, i')|,^,,Fk(t,t), (13) 

wp(i) = \Jdtdt'Fi,{t,t')l^^, I F^{t,t). (14) 

The evolution of the distribution n\^{t) and frequency 0Jp{t) is shown in Fig. ^ 
for the Hartree and Basketball approximations. In the Hartree case there is no 
inter- mode communication (no scattering), so the system does not equilibrate. In 
the Basketball case, the distribution and dispersion relation approach equilibrium 
profiles around mt w 80. To find out if the system has indeed reached a universal 
equilibrium state, we study the evolution of individual modes for different initial 
conditions with identical energy density (see Fig.|21). One can see that the individual 
modes become roughly equilibrated at mt « 1000. Further on there is practically 
no inter-mode communication, so it is reasonable to say that kinetic equilibration 
has been established. There is still, however, a slow drift towards the complete 
equilibrium state, for which the mode distributions should agree for the various 
initial conditions. As one can see in Fig.|21 the evolution from the initial conditions 
Tl and T3 differs substantially even at mt ~ 2000, so the system is still far from 
complete equilibrium. The reason for the slow drift towards full equilibrium at later 
stages seems to be the fact that the system is not chemically equilibrated yet. This 
can be seen by plotting the evolution of the total particle number density n^^i(t) = 
^k(i)j shown in Fig.|2| One observes that the kinetically pre-equilibrated state 
remembers the initial particle number, and chemical equilibration occurs at a very 
slow rate. 

Fitting the distribution nk(t) and frequency i^p{t) with an equilibrium form 
we can extract an effective temperature T„ff(t), chemical potential /ictf(^) and mass 
M^ff{t) (see Fig. 0)). Exponential fits for the effective temperature and chemical 
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Fig. 2. Evolution of individual modes Fig. 3. Evolution of total particle num- 
for various initial conditions with same ber density nt„t for (T1,T2 and T3) 
energy density (T1,T2 and T3). 
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Fig. 4. Evolution of effective mass, temperature and chemical potential (Tl case). 
Hartre mass Mh is included for comparison. 

potential indicate that final equilibrium is reached at mt ~ 10^~^. Comparing to 
the case in 2+1 dimensions [[7|, chemical equilibration seems to be much slower. 
The values of the effective masses seem to be roughly the same in the Basketball 
and Hartree approximations (also included in Fig. 0)). 

4. Broken Phase 

In the broken phase 0^0 and hence one can compare the approach to equilibrium 
between the two-loop and Basketball approximations. This is particularly interest- 
ing for the two-loop case, since perturbation theory predicts no on-shell scattering 
for the corresponding self-energy contribution. Thus a study of the time evolution 
using, for instance, Boltzmann-like equations will fail to show equilibration. An ap- 
proach based on the 2PI effective action, however, includes higher-order self-energy 
contributions that do contain on-shell scattering. It also takes into account off-shell 
scattering, which can play an important role far away from equilibrium. 





Kinetic and Chemical Equilibration in Scalar Theory 



7 




Fig. 5. Evolution in time of the occupation numbers, nk vs. Wk (left), and dispersion 
relation ui'^ vs. (right), for the Tl initial condition. We compare the results of 
the two-loop (green/gray) and the Basketball approximation (black). 
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Fig. 6. Evolution of total particle number density n^t for the two-loop and Bas- 
ketball approximations. 

The system is initialized with a distribution function of the form (|12|) . The 
mean field is taken at the tree level vacuum expectation value <j>(t = 0) = Wtj„e = 
a/6|to2|/A. For not so large couplings this is close to the actual variational value, 
so the time evolution of (/)(t) should not affect much the dynamics of the two-point 
functions that we use to monitor the process of equilibration. In our simulations 
we use A = 1.^ The evolution of the distribution function and dispersion relation 
defined according to (fU^ and l(Tl|l is shown in Fig. Surprisingly, the process 
of equilibration in the two-loop approximation seems to occur almost as fast as 
in the Basketball case. The same holds for the chemical equilibration that takes 
place at later stages of the evolution, as one can see in Fig. As in the symmetric 
case, chemical equilibration is very slow, which impedes the approach to complete 
equilibrium. 

^This requires a longer memory kernel. We cut the memory at mtcut = 84. 
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5. Conclusions 

In this work several truncations of the 2PI effective action have been used to study 
equiUbration in ip^ theory in 3+1 dimensions. Early in the evolution, the distri- 
bution function and dispersion relation stabilize and take equilibrium-like forms 
due to a relatively fast kinetic equilibration (see also [El)- This occurs for both 
the two-loop and Basketball approximations, which take scattering into account. 
Surprisingly, equilibration is equally rapid in both approximations. An analysis of 
damping rates of distrubances close to equilibrium leads to a similar conclusion [ 

cni. 

The pre-equilibrated state has not completely lost the memory of the initial 
state. In particular, initial conditions with different total particle number pre- 
equilibrate to different states. The later approach to complete equilibrium is very 
slow and mainly determined by the time scales associated with chemical equilibra- 
tion. 
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